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Abstract. The Japanese shareholding network existing at the end of March 2002 is studied empir- 
ically. The network is constructed from 2,303 listed companies and 53 non-listed financial institu- 
tions. We consider this network as a directed graph by drawing edges from shareholders to stock 
corporations. The lengths of the shareholder lists vary with the companies, and the most compre- 
hensive lists contain the top 30 shareholders. Consequently, the distribution of incoming edges has 
an upper bound, while that of outgoing edges has no bound. The distribution of outgoing degrees is 
well explained by the power law function with an exponential tail. The exponent in the power law 
range is y = 1 .7. To understand these features from the viewpoint of a company's growth, we con- 
sider the correlations between the outgoing degree and the company's age, profit, and total assets. 
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INTRODUCTION 

The economy is regarded as a set of activities of irrational agents in complex networks. 
However, many traditional studies in economics investigate the activities of rational and 
representative agents in simple networks, i.e., regular networks and random networks. 
To overcome the limitations of such an unrealistic situation, the viewpoint of irrational 
agents has emerged. However, simple networks have been adopted in economics and 
many studies of agent simulation. Recently^the study of complex networks has revealed 
the true structure of real- world networks ll|]|2|]. Gross networks such as interindustry 
relations have been considered in the macro economy, but the detailed structures of net- 
works constructed from entities in the micro economy have not been clearly elucidated 

By common practice, if we intend to discuss networks, we must define their nodes 
and edges. Edges represent the relationships between nodes. In this study, we consider 
companies as nodes. To define the relationships between companies, we use three view- 
points: ownership, governance, and activity. These three viewpoints have relations with 
each other. Ownership is characterized by the shareholding of companies, which is the 
subject of this article. Governance is characterized by the interlocking of directors, and 
it is frequently represented by a bipartite graph constructed from corporate boards and 
directors. The activity networks are characterized by many relationships: trade, collabo- 
ration, etc. 

In this article we consider Japanese shareholding network at the end of March 2002 
(see Ref. |@] for shareholding networks in MIB, NYSE, and NASDAQ). We use data 
published by TOYO KEIZAI INC. This data source provides lists of shareholders for 
2,765 companies listed on the stock market or the over-the-counter market. The lengths 



of the shareholder lists vary with the companies. The most comprehensive lists contain 
information on the top 30 shareholders. Types of shareholders include listed companies, 
non-listed financial institutions (commercial banks, trust banks, and insurance compa- 
nies), officers, and other individuals. In this article, we only consider the shareholding 
network constructed from 2,303 companies listed on the stock market and 53 non-listed 
financial institutions. Accordingly, the size of this network is N = 2,356, and the total 
number of edges is L = 22,435. 

This paper is organized as follows. In the next section we consider the degree distri- 
bution for both incoming edges and outgoing ones and show that the outgoing degree 
distribution follows a power law function with an exponential cutoff. In the following 
section, we discuss correlations between the degree and the company's age, profit, and 
total assets. This is because we assume that the dynamical change and growth of busi- 
ness networks can be explained by the company's growth. The last section is devoted to 
a summary and discussion. 

DEGREE DISTRIBUTION IN THE DIRECTED SHAREHOLDING 

NETWORK 

If we draw arrows from shareholders to stock corporations, we can represent a share- 
holding network as a directed graph. If we count the number of incoming edges and that 
of outgoing edges for each node, we can obtain the degree distribution for incoming de- 
gree, k[ n , and that of outgoing degree, fc ut- However, as explained above, the lengths of 
the shareholder lists vary with the companies, and the most comprehensive lists contain 
the top 30 shareholders. Therefore, the incoming degree has an upper bound, k m < 30, 
while the outgoing degree has no bound. 

The log-log plot of the degree distribution is shown in the left panel of Fig. [H and 
the semi-log plot is shown in the right panel of Fig. d In this figure, the horizontal 
axis corresponds to the degree and the vertical axis corresponds to the rank. The open 
circles represent the distribution of the incoming degree, and the plus symbols represent 
that of the outgoing degree. As mentioned above, we can find the upper bound for the 
incoming degree. In the discussion below, therefore, we consider only the distribution 
of the outgoing degree. 

The solid lines correspond to the fitting by the power law function p(k out ) k J t with 
the exponent y = 1.7, and the dashed lines correspond to the fitting by the exponential 
function. We can see that the outgoing degree distribution follows the power law distri- 
bution with an exponential cutoff. The exponential part of the distribution is constructed 
from 40 nodes, which are 38 financial institutions and 2 trading firms. On the other hand, 
the power law part of the distribution is constructed from 2,316 nodes, and 96% of them 
represent non-financial institutions. 

The above results suggest that different mechanisms work in each range of the dis- 
tribution. We assume that the dynamics in the range of the exponential distribution is 
essentially explained by the model proposed in Ref. ED. This model is based on a so- 
called BA model [8], i.e., the preferential attachment in growing networks that explains 
the power law distribution of the degree. Here, we assume that the outgoing degree dis- 
tribution in the tail part is explained based on the BA model, since financial institutions 
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FIGURE 1. Log-log plot (left) and semi-log plot (right) of the degree distribution . The open circles 
represent the distribution of the incoming degree, and the plus symbols represent that of the outgoing 
degree. The solid lines correspond to the fitting by the power law function with the exponent y = 1 .7, and 
the dashed lines correspond to the fitting by the exponential function. 



must invest money as shares. Therefore, financial institutions actively obtain shares of 
companies newly listed on the stock market or the over-the-counter market. This mech- 
anism is almost the same as preferential attachment. However, we must modify the BA 
model in order for it to explain the exponential distribution of the outgoing degree. If 
we consider an extended BA model with aging of the nodes, then we can obtain the 
exponential degree distribution [7]. However, as explained in the next section, this is not 
the case. On the other hand, if we consider an extended BA model that includes the cost 
of adding edges to the nodes or the limited capacity of a node, then we can obtain the 
exponential degree distribution 0]. As explained in the next section, this is actually the 
case. 

We consider cliques in networks to be important characteristics for constructing a 
model that can explain the power law distribution of the outgoing degree. Cliques in 
networks are quantified by a clustering coefficient 02D, which is defined in the case of 
undirected networks. Supposing that a node i has k[ edges, then at most — l)/2 
edges can exist between them. The clustering coefficient of node i, Q, is the fraction of 
these allowable edges that actually exist for e,-, i.e., Q = 2ej/ki(ki — 1). This has been 
calculated for an undirected shareholding network by using the same data used in this 



article [5], and it has been shown that the clustering coefficient follows the power law 
distribution, C(k) k~ a , with the exponent a = 1.1. Such a scaling property of the 
distribution of clustering coefficients has also been observed in biological networks, and 
this has motivated the concept of hierarchical networks 1 1 1 it 1 ill . We believe that the 
dynamics in the power law distribution of the outgoing degree is essentially explained 
by extending the model proposed in Ref. lillll . 

The clustering coefficient is approximately equal to the probability of finding triangles 
in the network. The triangle forms the minimum loop. Therefore, if node i has a small 
value of Cj, then the probability of finding loops around this node is low. Consequently, 
this scaling property of the clustering coefficient suggests that the network is locally 
tree-like. 

CORRELATION BETWEEN OUTGOING DEGREE AND 
COMPANY'S GROWTH 

We consider the correlations between the outgoing degree and the company's age, profit, 
and total assets. We believe that knowing the characteristics of nodes is useful for 
constructing models explaining the dynamics of networks. In many complex networks, 
it is difficult to quantitatively characterize the nature of nodes. However, in the case of 
economic networks, especially networks constructed from companies, we can obtain 
the nature of nodes quantitatively based on balance sheets, and income statements, 
for example. We consider this a remarkable characteristic of business networks, and it 
allows us to understand business networks in terms of the company's growth. We believe 
that the dynamics of business networks must be explained by the theory of company 
growth. 



Outgoing degree and company's age 

As mentioned in the previous section, we consider the correlation between the outgo- 
ing degree and the company's age in order to clarify whether the BA model with aging of 
the nodes is applicable to a particular case. In the following, we measure the company's 
age in units of months. 

The log-log plot of the distribution of the company's age is shown in the upper left 
panel of Fig. El In this figure, the horizontal axis is the age in the unit of months, and 
the vertical axis is the rank. This figure shows that the distribution does not fit a simple 
function such as a power law function. 

The semi-log plot of the distribution of the company's age is shown in the lower left 
panel of Fig. [2j In this figure, the meaning of the axes is the same as in the upper figure. 
The dashed line corresponds to the fitting by the exponential function. This means that 
the distribution of the company's age approximately follows the exponential distribution. 
It is expected that the age of a company has a relation with its lifetime, and the lifetime 
of bankrupted companies follows an exponential distribution lfl2lllfl3ll . 

The upper right panel of Fig. El shows a log-log plot of the correlation between the 
company's age in months and the outgoing degree that follows the exponential distri- 
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FIGURE 2. Log-log plot (upper left) and a semi-log plot (lower left) of the distribution of the company's 
age. The dashed line corresponds to the fitting by the exponential function. Log-log plots of the correlation 
between the company's age and the outgoing degree that follows the exponential distribution (upper right), 
and the one that follows the power law distribution (lower right). 



bution. We can observe that these two quantities have no correlation. To quantify this 
observation, we calculated Spearman's rank correlation S and Kendall's rank correlation 
K and obtained S = 0.27 and K = 0.17, respectively. These results mean that there is 
no correlation between the company's age and the outgoing degree that follows the ex- 
ponential distribution. This suggests that the BA model with aging of the nodes is not 
applicable to this case. 

The lower right panel of Fig. shows a log-log plot of the correlation between 
the company's age in months and the outgoing degree that follows the power law 
distribution. We can also observe that these two quantities have no correlation. In this 
case, S = 0.32 and K = 0.24. These results mean that there is no correlation between the 
company's age and the outgoing degree that follows the power law distribution. Needless 
to say, the BA model with aging of the nodes is not applicable to this case. 

A comparison of these two figures clarifies the independence between the outgoing 
degree and the company's age. There are some companies with a few outgoing edges in 
the range where the lifetime is longer than 10 3 months. 



Outgoing degree and company profit 



We consider the correlation between the outgoing degree and the company's profit. 
Here the company's profit is the amount of money that the company gains when it is 
paid more for something than it cost to make, get or do it. Hence the profit includes 
information on the flow in the network. Here, we consider the company's profit as an 
amount of money stored in the period from the beginning of April 2001 to the end of 
March 2002. 

The log-log plot of the distribution of the company's profit is shown in the upper left 
panel of Fig.|3j In this figure, the horizontal axis is profit / in units of million yen, and 
the vertical axis is the rank. The solid line corresponds to the fitting by the power law 
function. This is represented by the probability density function (pdf), p(I) as p(I) I a 
with the exponent a = 2. This shows that the distribution in the profit in the middle 
range follows the power law distribution. The dashed line corresponds to the fitting by 
the exponential function. The exponential distribution in the high-profit range is better 
clarified by the semi-log plot. 

The semi-log plot of the distribution of the company's profit is shown in the lower left 
panel of Fig.|3l The meanings of axes, solid line, and dashed line are the same as those 
in the upper panel. This figure shows that the profit in the tail part follows an exponential 
distribution. However, the number of companies within this range is small. 

The upper right panel of Fig. |3] shows a log-log plot of the correlation between the 
company's profit in the unit of million yen and the outgoing degree that follows the 
exponential distribution. As mentioned previously, the exponential part of the outgoing 
degree distribution contains 40 nodes, which are 38 financial institutions and 2 trading 
firms. Almost all of these financial institutions are not listed on the stock market, and we 
could not obtain their profit data. Needless to say, there is also a possibility that the profit 
of these financial institutions was negative. Therefore, the total number of data point is 
less than 40. We can see that these two quantities have no or only weak correlation. 
Spearman's rank correlation and Kendall's rank correlation are S = 0.31 and K = 0.23, 
respectively. These results mean that there is no or only weak correlation between the 
company's profit and the outgoing degree that follows the exponential distribution. 

The lower right panel of Fig. |3] shows a log-log plot of the correlation between the 
company's profit and the outgoing degree that follows the power law distribution. We can 
see that these two quantities also have no or only weak correlation. In this case, S = 0.32 
and K = 0.24. These results mean that there is no or only weak correlation between the 
company's profit and the outgoing degree that follows the power law distribution. 

Here, we consider the profits of the companies comprising the network. Therefore, the 
distribution of profit, i.e., the upper left panel of Fig. O is not so impressive. However, 
the profit of Japanese companies has remarkable properties: (i) It is confirmed that in 
the period after 1970 the high-profit range followed the power law p(I) / , i.e., 
Pareto-Zipf's law, and was stable ifbUl . (ii) It is confirmed for the years 2000 and 2001 
that detailed balance was maintained for the top 70,000 companies in the high-profit 
range I15I1 . (iii) It is also confirmed for the years 2000 and 2001 that Gibrat's law was 
maintained for the top 70,000 companies in the high-profit range fl5ll . Here, Gibrat's 
law means that the growth rate of profit is independent of the profit in the initial year, 
(iv) These characteristics are related to each other I15I1 . 
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FIGURE 3. Log-log plot (upper left) and a semi-log plot (lower left) of the distribution of the company's 
profit. The solid lines correspond to the fitting by the power law distribution with the exponent a = 2, and 
the dashed lines correspond to the fitting by the exponential function. Log-log plots of the correlation 
between the company's profit and the outgoing degree that follows the exponential distribution (upper 
right) and the one that follows the power law distribution (lower right). 



Outgoing degree and company's total assets 

We next considered the correlation between the outgoing degree and the company's 
total assets. The company's total assets are naively regarded as the amount of profit. A 
portion of the company's total assets is invested in shares. Accordingly, it is expected 
that the company's total assets have a relation with the cost of adding edges to the nodes 
or with the limited capacity of a node. 

The log-log plot of the distribution of the company's total assets at the end of March 
2002 is shown in the upper left panel of Fig. |4] In this figure, the horizontal axis is the 
total assets A in the units of million yen, and the vertical axis is the rank. The solid 
line corresponds to the fitting by the power law function p(A) A~ b with the exponent 
b = 1.6. This figure shows that the level of total assets in the middle range follows the 
power law distribution. The dashed line corresponds to the fitting by the exponential 
function. 

The semi-log plot of the distribution of the company's total assets is shown in the 
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FIGURE 4. Log-log plot (upper left) and a semi-log plot (lower left) of the distribution of the company's 
total assets . The solid lines correspond to the fitting by the power law distribution with the exponent 1 .6, 
and the dashed lines correspond to the fitting by the exponential function. Log-log plots of the correlation 
between the company's total assets and the outgoing degree that follows the exponential distribution 
(upper right), and the one that follows the power law distribution (lower right). 



lower left panel of Fig.|4| The meanings of axes, solid line, and dashed line are the same 
as in the upper panel. This figure shows that the level of the total assets in the tail part 
of the distribution follows the exponential distribution. However, the change from the 
power law distribution to the exponential distribution is not smooth. 

The upper right panel of Fig. HI shows a log-log plot of the correlation between the 
company's total assets at the end of March 2002 in units of million yen and the outgoing 
degree that follows the exponential distribution. Contrary to the case of the company's 
profit, we can obtain data on the company's total assets for the non-listed financial 
institutions. This figure shows that these two quantities strongly correlate. Spearman's 
rank correlation and Kendall's rank correlation are S = 0.71 and K = 0.55, respectively. 
These results mean that there is strong correlation between the company's total assets 
and the outgoing degree that follows the exponential distribution. This suggests that the 
BA model incorporating the cost of adding edges to the nodes or the limited capacity of 
a node is applicable to this case. 

The lower right panel of Fig. |4| shows a log-log plot of the correlation between the 



company's total assets and the outgoing degree that follows the power law distribution. 
We can see that these two quantities also have strong correlation. In this case, S = 0.65 
and K = 0.5l. These results mean that there is strong correlation between the company's 
total assets and the outgoing degree that follows the power law distribution. 

This suggests that the level of the company's total assets has strong correlation with 
the outgoing degree. In particular, it is assumed that the outgoing degree distribution in 
the tail part is explained by the BA model incorporating the cost of adding edges to the 
nodes or the limited capacity of a node. However to confirm this assumption, we must 
study the accumulation process of the company's total assets. 

SUMMARY 

In this paper we considered the Japanese shareholding network existing at the end of 
March 2002. Although we could not clarify the detailed characteristics of the incoming 
degree distribution, we could obtain useful information on the outgoing degree distri- 
bution: (i) The outgoing degree distribution follows the power law with an exponential 
cutoff, (ii) The important factor in the growth of the business network is not the com- 
pany's age but its total assets. This means that old companies do not necessarily have 
large total assets, (iii) It is expected that the dynamics of the tail part of the outgoing 
degree distribution, i.e. the exponential distribution range, can be explained by the BA 
model incorporating the cost of adding edges to the nodes or the limited capacity of a 
node. However, the mechanism for the power law part of the outgoing degree distribu- 
tion is still not known. We believe that the dynamical change and growth of business 
networks must be explained by the company's growth. Consequently, knowing the dy- 
namics of the company's growth is a key concept in considering the growth of economic 
networks 00. 

We would like to conclude with two observations. The first concerns degree correla- 
tion. In this article, we considered the shareholding network as a directed graph. How- 
ever, if we ignore the direction of edge, we can calculate many quantities to obtain the 
characteristics of networks. For example, the nearest neighbors' average degree of nodes 
with degree k, (k nn ) (k), is an important quantity W18I1 . This was calculated by Ref. \3\ for 
the Japanese shareholding networks at the end of March 2002, in which (k rm )(k) k~ v 
was obtained with the exponent v = 0.8. This means that hubs are not directly connected 
to each other in this network, i.e., it's a degree non-assortative network. 

The second observation concerns the spectrum of the graph. It has recently been 
shown through an effective medium approximation that the probability density function 
of the eigenvalue for the adjacency matrix, p(A), is asymptotically represented by that 
of the degree distribution p(k), i.e., p(A) ~ 2\X\p(X 2 ), if the network has a local 
tree-like structure lfl9ll . As mentioned previously, the shareholding network has this 
characteristic. Therefore, if p(k) asymptotically follows the power law distribution, 
p(A) also asymptotically follows the power law distribution. In addition, we can derive 
the scaling relation 8 = 2y — 1 . At the end of March 2002, the tail part of the eigenvalue 
for the adjacency matrix of the undirected Japanese shareholding network followed 
p(A) with the exponent 8 = 2.6. On the other hand, the degree distribution 

in this case is p(k) k~? with the exponent y = 1.8. Therefore, the scaling relation is 



guaranteed 
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